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Abstract: The chemical potential of the electron gas on a two-dimensional 
rectangular lattice is determined. An approximate expression for exp(— n/T) 
is obtained, and its second order approximation is discussed in some detail. 
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Introduction 



Investigations of correlated electron systems of various dimensionality are 
one of the priorities of contemporary condensed matter physics (for example 
[l]).The driving force behind the scientific interest in such systems is the 
fact that they exist in high temperature superconductors and organic con- 
ductors. Although the existence of correlated electrons in these materials is 
regarded as an "established fact", it is far from being completely clear how- 
does one theoretically describe them [2]. The question is if it is possible to 
use the Fermi liquid model, or the Luttinger liquid model has to be used. For 
a recent introduction to the Luttinger liquid theory see for example [3]. 

The aim of the present paper is to calculate the chemical potential of the 
electron gas on a two dimensional rectangular lattice. Theoreticallythis is an 
interesting problem in itself,and the experimental motivation is amply dis- 
cussed in the literature (for example [4]). The calculation will be performed 
within the standard Fermi liquid (FL) model. This presumption could be crit- 
icized on the grounds that there are recent results indicating the breakdown 
of the FL model in cuprate superconductors and heavy-fermion systems(such 
as [5]). The idea of the present calculation is to contribute to this result, in 
a different way: by calculating the chemical potential,then separately the 
electrical conductivity and finally drawing conclusions by comparison with 
experimental data. 

Calculations 

The lattice axes are designated by x and y.It can be shown that the number 
density of a 2D Fermi liquid is given by [6]: 

_ <N > r exp{(3{n - e)) d 2 p 
U S J 1 + exp(/3((j, - e)) (2nhf { ' 

where n denotes the number density and all the other symbols have their 
usual meanings and p = hk. It then follows that dp = h^de, where e denotes 
the energy.Accordingly, 

d p = dp x dp y = h — —^-de x de y (2) 
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Assuming that the energy bands along each of the lattice axes have the 
form e = —2tcos{ks), where t denotes the hopping integral and s the lattice 
constant along the axis, it follows that 



dk x 1 de q 



de x 2s x t x 7(1 - (e x /(2t x )f) 
An analogous relation is valid along the y axis. This implies that 



(3) 



(f — ^ de x de y 

4s x t x s y t v v/(l - (e x /(2t x )) VI - (e,/(2t s/ )) 2 1 ; 

Using the expressions for the particle energy along the lattice axes,the expo- 
nential factor in eq.(l) can be expressed as 

exp(P(fi - e)) = exp(/3(/i - (e x + e y ))) (5) 

Inserting eqs.(4) and (5) into eq.(l), after some algebra leads to 



n 



8tt 2 s x t x s y ty J J 1 + exp((3(e x + e y - //)) 
1 1 
7(1 - (e x /(2t x )*) 7(1 " (^/(2t,) 2 )) 



(6) 



which is the required equation linking various measurable parameters of 
the electron gas and the lattice with the chemical potential of the electron 
gas. The limits of integration are ±2t along each of the lattice axes, and the 
problem is how to solve this integral, which can not be done in closed ana- 
lytical form. Looking at things in a purely mathematical waythe easiest way 
to solve this integral is to develop the sub-integral function into series. Rea- 
soning in this waythe Fermi distribution function can be re-expressed as 
follows, under the condition that e x + e y -< /i: 

1 l = oo 

Introducing this development into eq.(6),and slightly grouping the terms 
in it, one gets the following expression: 
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■y 1=00 ,. ,. 

n = s T 2 c ; " ; 12 / / expW^x + £y- A*)) 

OTT O X i X 0yiy J J 

de x de y 

7(1 - (e^k) V(l " (^/2t„) 2 ) ( } 

Integrating this expression, within the limits of ±2t along both axes, and 
performing some algebra, leads to the following final result for the filling factor 
of a Fermi gas on a two dimensional lattice: 

1 oo 

n = ^(-l) i exp(-/3/i/)/ (2/t^)/o(2/^) (9) 

where the symbol I n (x) denotes the modified Bessel function of the first 
kind of the order n and (5 is the inverse temperature. The functions I n (x) are 
defined by: 

fe=oo ( /n\n+2k 

Discussion 

Equation (9) represents the general result of this paper, expressed in its 
most compact form. The obvious question is the applicability of this result to 
real physical systems. Limiting the summation in eq(9) to terms with / < 2 
gives 

n = ^-[l-exp(- / ,/T)/ (^))/ ^ + 

exp(-2 / i/T)/ (^)/o(^)] (11) 
Solving this equation for exp(— fi/T) leads to 
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exp(-n/T) = (/ (^)/o(^) ± 



^(f^fO ~ 4(1 - 2n S , S ,)/ (^)/ (|^))/(2/ (^)/ (^)) (12) 

Developing this equation into series up to second order terms in t x and 
t y as small parameters, gives the following expression for exp(—/i/T): 

exp(-f,/T) = i(l - y^s, - 3) + \[{t x /Tf + (t,/T) 2 ] 
, 16^5^ - 5 _ 3 ^ + y^^j + 16^^(9-167^3,,) - 17 ^ + ^ 



v y/Sns x Sy - 3 7 2 V T 2/V (8ns x s„ - 3) 3 / 2 

Compare this result with [7] where an expression for the chemical poten- 
tial of a ID electron gas has been derived.lt follows from eq.(13) that there 
exists a lower bound on the values of n for which it is applicable. The con- 
dition for validity of eq.(13) is 8ns x s v — 3^0 which implies n >- g 3 .The 
implication of this limitation is that eq.(13) can not be used for extremely 
small values of n. Inspection of this expression shows that it represents the 
dependence of the chemical potential of a 2D electron gas on various mea- 
surable parameters of the system. As the lattice constants along both lattice 
axes are taken into account, it follows that this equation gives the possibility 
of analyzing the influence of high external pressure on the chemical potential. 
This possibility could turn out as being useful in calculations of the pressure 
dependence of the electrical conductivity of 2D organic conductors. 

An interesting problem is the value of n (the filling factor) for which the 
chemical potential becomes equal to zero. It follows after some algebra that 
fj, = for n = no, where 

^ 7T 4_ T2[1 + 16(t 2 +t 2 )]+32(t 2 +t 2 )2 

l Qs x s y {T^-2(tl + tl)Y 

jT* + 10T*{tl + + & + iff (14) 
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The outright implication of this result is that the Lieb-Wu theorem [8] 
is not generally applicable to 2D systems. It can be applied only for certain 
values of various parameters of the system which enter into eq.(14).It also 
follows from the last equation that no is a function of the external pressure 
to which the specimen is subdued. 

Conclusions 

In this paper we have determined the chemical potential of the electron 
gas on a 2D rectangular lattice. A general expression linking the filling fac- 
tor, the chemical potential and various measurable parameters of this system 
was obtained. Starting from this general expression, an approximate expres- 
sion for the chemical potential has been obtained. It turns out that this 
expression is much more complicated than the corresponding expression for 
a ID electron gas. This expression can be applied in analyzing real experi- 
mental data, and in calculations of the electrical conductivity of Q2D organic 
metals,which will be done in future work. 
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